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As it was proposed and recently verified experimentally, the mechanism of Feshbach resonance in 
a condensate can create a second condensate component of molecules that coexists with the atomic 
condensate. In this work we investigate signatures of the presence of the molecular condensate 
through the equilibrium properties and collective excitations of the hybrid system of atoms and 
molecules, subjected to a trap, employing a time-dependent variational ansatz. We show that the 
shape of the condensate changes significantly by the presence of the molecules and that modes 
unique to this hybrid system appear at observable frequencies. 



The recent realization of Bose-Einstein condensation 
(BEC) in a dilute gas of alkali atoms has open a new 
opportunity for the theoretical and experimental investi- 
gation in quantum degenerate many-body systems Q, • 
In contrast to superfluid helium, these weakly interacting 
gases are much more amenable to theoretical prediction 
and quantitative analysis. After the observation of such 
condensates, a new generation of more complex exper- 
iments has been done involving the production of mul- 
tiple species BEC. Such mixture may consist of differ- 
ent hyperfine states of the same atom or different spin 
state of the same atomic species. In the former case 
the composition of mixture remains fixed and in the 
latter case the internal conversion leads to important 
new phenomena. It has been proposed that the interac- 
tion through the Feshbach resonance that bring a binary 
atom system to an intermediate state molecule recently 
observed in BEC H|, create in the dilute atomic Bose- 
Einstein condensate (BEC) a second molecular conden- 
sate component 0, @. We should mention also the case 
of a dilute atomic BEC coupled to a molecular Bose gas 
by coherent Raman transitions, which may lead to what 
is coined superchemistry @,IU. 

Recently, several groups have been working in order 
to detect a signature of molecule formation in BEC with 
Feshbach resonance Since these systems allow 

experimentalists to change the magnetic filed (and con- 
sequently the detuning) one can probe these systems in 
two different ways: 1 - adiabatic, where the detuning 
can be changed slowly (waiting for the system to equi- 
librate) and 2 - non-adiabatic, where a sudden change 
of the magnetic field is performed (usually with a pulse). 
The non-adiabatic change leaving the system out of equi- 
librium will produce oscillations, among other effects, in 
the number of atoms, and that was indeed observed 

In the adiabatic change the system is always at the 
equilibrium configuration. In this paper we show that in 
this regime the shape of the equilibrium condensate cloud 
of atoms and the behavior of the excitation energies of 



the collective modes can be a signature of the presence 
of molecules. It will be shown that the effects of the 
presence of molecules are more important and could be 
detected experimentally for weak resonances 

We start by reviewing the many body theory for a 
condensate on a Feshbach resonance. In a two component 
picture [f| the mean field hamiltonian density describing 
the trapped system is 
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°a(m) 



denote the macroscopic wave function of 



the atom (molecule) condensate, e is the detuning pa- 
rameter, a is the Feshbach resonance parameter in- 
troduced in Ref.[f|, and we assume that the confin- 
ing potential for the atoms and molecules are harmonic 
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with axial symmetry, V a = ^^{P^x 2 + /3^y 2 + /?"z 2 ), 
V m = muj^iP^x 2 + [3™y 2 + (3™z 2 ). As usual, the Gross- 
Pitaevskii equations of motion for this model can be de- 
rived by imposing that the action 
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where E[ip a ,tp m ] = J d 3 r Ti. , is stationary with respect 
to the variation of of <fi a ( m ), subject to the normaliza- 
tion condition constraint j d 3 r(\<j> a \ 2 + 2\<j> m \ 2 ) = N. By 
looking for solutions of the form (f) a (r,t) — e~ li z~ (f) a (r) 

and (j) m (r,t) = e~ l ~%~ (f> m (r) we find the coupled time- 
independent GP equations, 

h 2 

^ a = (-^V 2 + K + A Q |0 a | 2 + A am |</> m | 2 )(/> a 
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Now we turn our attention to the limit where the two- 
component theory reduces to an effective one-component 
theory. This limit is reached when the detuning is large 
enough that e greatly exceeds the kinetic energy as well as 
any interaction energies, i.e., e >> X m n m , X am n a , X a n a , 
where n a and n m are the atomic and molecular densi- 
ties respectively,plus the additional hypothesis that the 
number of molecules are small compared to the number 
of atoms n m /n a << 1. These assumptions allow us to 
neglect the interaction terms and the chemical potential 
(fi rs X a n a ) in the Gross-Pitaevskii equation 10}, leading 
to a simple relation between the molecular and atomic 
condensate 
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The insertion of this expression into equation © yields 
an effective single condensate Gross-Pitaevskii equation 
for 6 a 



A*0a = 
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Note that, from JSJ, the limit |</> m |/|</> a | << 1 implies 
e >> a(f> a - As A a is positive and we want to vary e 
in a way that X e f t is always positive, it follows from 
(Q that we must set e > a 2 /X a . Assuming that all the 
interaction strengths are of the same order of magnitude, 
the conditions of validity of the one-component theory 
can be summarized by the inequalities: e » X a n a ,e » 
OL\fna with e > a 2 / X a - In what follows we will consider 
two regimes characterized by the parameter r, 
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which relates the characteristic energies of our system, 

a 2 /X a =ra^/n a \,a^/n a ~ = rX a n a (9) 

i) Strong resonance limit- When r » 1 the character- 
istics energies satisfy the set of inequalities, a 2 /X a >> 
cty/n^ » X a n a and for the entire range of variation of 
the detuning, even up to A e // rj 0, the conditions of va- 
lidity of the effective one-component theory is satisfied. 
Therefore, when the parameter r satisfies the above in- 
equality one can picture the system as being one conden- 
sate of atoms with an effective scattering length. We will 
show, in our calculations that with values of r in this 



range, the two component and the effective one com- 
ponent theory agree even when the effective scattering 
length is almost zero, the molecules are practically inex- 
istent in this regime. 

ii) Weak resonance limit -On the other hand when r w 
1 it follows that all the characteristics energies are of the 
same order of magnitude a 2 / X a ~ a^/n^ ~ X a n a and for 
e w a 2 /X a ,the conditions of validity of the effective one- 
component theory is not satisfied. This is corraborated by 
our calculations where we show that, for detunings of the 
order of a 2 /X a that correspond, in the single condensate 
effective theory, to almost zero effective scattering length, 
we obtain significantly different equilibrium configuration 
and also have unique modes to the hybrid system that 
are relatively low in energy (u> e xc ~ ^>trap)- 

The weak resonance regime proves to be an ideal place 
to experimentally detect a clear signature of molecule 
formation. The expression for the parameter r illustrates 
that there are three ways of approaching this regime: by 
increasing the condensate density, having systems with 
a higher asymptotic atom-atom scattering length and/or 
reducing the value of a. 

Next we will show, as pointed out above, that the equi- 
librium spatial distribution of the atoms and the excita- 
tion energies of specific collective modes have different 
behaviors as a function of the detuning e depending on 
the regimes discussed above and can be a signature of 
the presence of molecules in the system. 

The standard approach to investigate the properties 
of a hybrid system would be to solve the equations J3J), 
(0} and © to find the equilibrium shapes and to solve 
the linear approximation of its time dependent extension 
to find the collective excitations. However in this paper 
we use a simpler approach, the variational method. In 
this method we employ a Gaussian ansatz 0, to pa- 
rameterize the time dependence of the condensate wave 
function as 

cj> a{m) (x, t) = A a{m) (x, t)e lF ^> (10) 

where the amplitudes of the condensate wave functions 
are given by 



A a = 
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and the phase as 
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.3=1,2,3 ^ 

The variational parameters qj {Qj), r| (Rj) correspond 
respectively to the widths and translation of the atomic 
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and molecular condensate clouds, which have an ellip- 
soidal shape pj(Pj), TTj (IE) are the corresponding mo- 
menta. The N a ( m ) are the number of particles of each 
component 9 a / m \ the corresponding phase. 

When we replace the condensate wave function 
parametrized as in eas. (|1()I14|I into equation (J2J), T re- 
duces to a classical action 



a -PjQj 



., = 1.2.3 
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where the energy E has the form E = E a +E m +Ei+E FR , 
with Evj( m ) being the usual Gaussian Gross-Pitaevskii en- 
ergy for a single atomic(molecular) condensate 0, Ei 
are the interaction energies between these two species 
and Efr is the Feshbach component. 

The variation of the action with respect to the twenty 
eight parameters leads to hamiltonian type equations of 
motion for the time evolution of these parameters. De- 
noting collectively the parameters by Xi, i — 1, • • • , 28, 
the equations of motion can be written as 01 



£r w (x)x, 



dE 
dX~k 
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where (X) is the anti-symmetric matrix coming from 
the variation of the "kinematic" piece of the action (14) 
E(X) depends only on the phase difference, 9 = 9 a — 
9 m /2 and as a consequence the quantity TV = N a + 2N m 
is a constant of the motion. This reduces the number 
of degrees of freedom to thirteen which we take as the 
relative phase 9 = 9 a — 9 m /2 and the relative population 
n = (N a — 2N m )/2, besides the translational and shape 
degrees of freedom. 

i) Equilibrium configuration 

The equilibrium configuration is determined by the 
condition = 0, where k = 1, 26 which leads to 



dE 
~dX~k 



(X ) - , 



(17) 



where Xo represents the equilibrium values of the pa- 
rameters. Thirteen of the equations l|17H lead immedi- 
ately to pj = 0{Pj = 0), tt/ = 0(IEJ = and 9 = tt/2 
which is the phase of the minimum configuration. Six of 
the remaining thirteen equations refers to the equilibrium 
displacement of the condensate clouds. These equations 
have the trivial solution r| = 0(Rj = 0),j = 1,2,3, which 
is the one considered in this paper. Therefore we are left 
with seven equations to find the equilibrium widths of 
the condensate clouds and of the relative population n, 
n = (N a - 2N m )/2. 

We proceed by presenting the results of the static cal- 
culations obtained by solving (|17|) in the weak and strong 



FIG. 1: Plot showing the atomic condensate width as 
a function of the detuning in the weak resonance limit 
(r ~ 1) for (a) the radial width q x = q y = q r (b) the axial 
width q z . The dashed lined corresponds to the single com- 
ponent effective theory and the continuous line the two 
component molecular theory. The values of parameters 
used are: u) a = = 2tv x 500Hz, f3 r — 1, f3 z 
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resonance limits respectively. In the strong resonance 
limit we used r ~ 10. We computed the equilibrium ra- 
dial and axial widths (q x , q y , q z ) for the atomic conden- 
sate for different values of the detuning. We varied the 
detuning from 27rxl05 MHz to 27rx22 MHz correspond- 
ing, in the effective theory, to a A e / / varying from 80% to 
5% of the asymptotic interaction strength, respectively. 
The results are shown on Table I. 

One can see no difference between the two component 
molecular theory and the single condensate effective the- 
ory. This shows that all attempts, through an adiabatic 
change of the detuning, to detect molecular signatures 
in this regime will fail, because for all practical pur- 
poses the number of molecules formed is quite insignifi- 
cant. To achieve the weak resonant limit we used r ~ 1. 
The detuning values were changed from 2ttx41.7 KHz 
to 27rx8.7 KHz corresponding to the same effective in- 
teraction strength as in the strong resonance case (from 
80% to 5% of the asymptotic interaction strength). The 
results, shown on Table II, indicate a very different be- 
havior between the molecular and the single condensate 
effective theories, when the detuning approaches a 2 /\ a - 

In the two component theory, the size of the atomic 
condensate will not change significantly when we vary 
the detuning, due to the presence of molecules in contrast 
with a drastic change in the single condensate effective 
picture. We conclude that this regime proofs to be ideal 
for probing molecular signatures provided the detuning is 
adiabatically changed. The weak resonance limit results 
are displayed in Fig.l. 
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ii) Excitations 

The excitations are associated with the normal modes 
of the linear approximation of the equations of motion 
1(15)1. In this limit (jTp|l reduces to [3 

^2r kl (x )m =J2 n M( x o)m . (is) 
i i 

where Tlki(^-o) is the Hessian matrix at equilibrium 

Kki — Qx k dXi (X°) an< ^ are the displacement from 
the equilibrium value Xi = Xiq + rji. The normal 
modes correspond to oscillatory solutions of l|18|) . and 
the eigenfrequencies are given by the characteristic equa- 
tion det(irSl — H) =0. The roots of the caracteristic 
equation are real and come out in pairs (O, — Q), each 
pair associated to one normal mode. 

In the single-condensate case we have six normal modes 
which are related to the translations of the condensate 
clouds and to shape oscillations [l^- For the later, two of 
these correspond to axially symmetric oscillations where 
the axial and radial oscillation are in phase (mode c) and 
the other where they are out of phase (mode b), the third 
one being non-axially symmetric with the oscillations in 
x and y directions out of phase and the axial direction at 
its equilibrium value (mode a). 

In the two component case the normal modes will be 
linear combinations of the modes described above for 
each component plus one mode related to the dynam- 
ics of population exchange between atoms and molecules, 
giving a total of thirteen normal modes. In what follows 
we will investigate how the eigenfrequencies depend on 
the Feshbach resonance parameter a and the detuning e. 

In our calculation we obtain the following general re- 
sults, independent of the value of r: a)The six transla- 
tional degrees of freedom are decoupled from the shape 
and number degrees of freedom, originating three atomic 
and molecular translations in phase (center of mass 
modes) and three out of phase (dipole modes). b)The 
remaining seven modes can be divided into two groups. 
The hrst group, where the shape and number oscillations 
are decoupled, consists of two modes which correspond 
to an in phase and an out of phase atomic and molec- 
ular oscillations where both components are in mode a. 
The second group consists of coupled shape and number 
oscillations with each component in modes b and/or c. 
c) In the single component theory we have six modes, 
three translational and three shape oscillations. The 
three translational modes correspond to the three cen- 
ter of mass modes of the hybrid system and the three 
shape oscillations correspond to the three lowest energy 
modes not counting the center of mass modes. 

For the strong resonance limit, r « 10, the three lowest 
energy modes, disrecarding the translational ones, for the 
hybrid system correspond to in phase atomic and molec- 
ular shape oscillations (b+b), (a+a) and (c+c) respec- 
tively, and have the same behavior in both theories as 
exemplified on Table I for the mode (a+a). 

In the weak resonance limit r w 1 of the three modes 
only the mode (a+a) behaves differently depending on 



FIG. 2: Plot comparing two of the mode energies, in 
units of the radial trap frequency, as a function of the 
detuning. The dashed lines correspond to the single com- 
ponent effective theory and the continuous line the two 
component molecular theory. The lowest mode and four 
highest modes are not shown in the figure. See caption of 
figure 1 for the parameters. 

2.2 | — , — , — , — ,— , — ,— , — , — , — , — , — ^ — ^ — , — , — ^ — ^ — I 

2.1 - 



2 




1.7 - \ 



1.6 - \ 

1.5 - — 

1.4 ' — ' — ' — ' — — ' — — — ' — ' — 1 — ' — — — — ' — ' — ' — — — ' 
8 13 18 23 28 

e (2ji KHz) 



which theory is used. In the single component theory 
the a mode approaches abruptly to the c mode when 
A e // — > 0, because of the degenerescence present in the 
ideal gas limit. On the other hand for the hybrid the- 
ory this abrupt change is not present. The absence of 
a near degenerescence is a signature of the presence of 
the molecules. The values of this frequency are shown on 
Table II and two of the three lowest energy modes, for 
both single and hybrid theories are displayed in Fig. 2. 

Another possible signature of a molecular condensate 
is the presence of the dipole mode frequencies in the two 
componet theory. The frequency of the center of mass os- 
cillation motion of the hybrid BEC is entirely determined 
by the frequencies of the trap and does not depend on 
the interatomic interaction Q. On the other hand the 
frequencies of the dipole oscillation of the hybrid system 
depend on the Feshbach resonance interaction parameter, 
a, and on the detuning parameter e. Since this mode is 
unique to the hybrid system its presence is a clear sig- 
nature of the presence of molecules. This being an in- 
teresting phenomenon to be explored experimentally. Of 
the three dipole modes the radial ones are degenerate 
and have a higher frequency than the axial one. Using 
the same parameters for the weak limit we plot in Fig. 3 
the frequency of this mode as a function of the detun- 
ing. From this results we see that the frequency of this 
mode, unique to the hybrid system, decreases when the 
detuning approaches a 2 /X a reaching observable values. 

In conclusion, we have shown, that provided we work 
in the weak resonance limit, the hybrid atom molecule 
BEC system behaves qualitatively differently from a sin- 
gle condensate with an effective scattering length, ex- 
hibiting different behaviours for the equilibrium shapes 
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FIG. 3: Plot showing the dipole frequency, in units of the 
radial trap frequency, as a function of the detuning for the 
hybrid condensate. The lower curve corresponds to an 
axial oscillation and the upper curve radial oscillations. 
The parameters are the same as in figure 1. See text for 
details 
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TABLE I: Comparison between the single condensate theory 
and the hybrid theory in the strong resonant limit. The Fes- 
hbach resonance parameter in this case is a = 510 4 /im 3//2 Hz 
and the other parameters are the same as in figure 1. 



and collective excitations as a function of the detuning. 
Besides for reasonable values of the physical parameters 
it appears a unique mode of the hybrid system whose 
frequencies are found to be of the order of the trap fre- 
quencies for detuning values of the order of KHz making 
the observation of this mode experimentally feasible. 
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